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Domain Switching Kinetics in Disordered Ferroelectric Thin Films 
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We investigated domain kinetics by measuring the polarization switching behaviors of polycrys- 
talline Pb(Zr,Ti)03 films, which are widely used in ferroelectric memory devices. Their switching 
behaviors at various electric fields and temperatures could be explained by assuming the Lorentzian 
distribution of domain switching times. We viewed the switching process under an electric field as 
a motion of the ferroelectric domain through a random medium, and we showed that the local field 
variation due to dipole defects at domain pinning sites could explain the intriguing distribution. 

PACS numbers: 77.80.Fm, 77.80.Dj, 77.84.Dy 



Domain switching kinetics in ferroelectrics (FEs) un- 
der an external electric field E ext have been extensively 
investigated for several decades jl|, d i, i, i, @, 0, 0, H|. 
The traditional approach to explain the FE switching 
kinetics, often called the Kolmogorov-Avrami-Ishibashi 
(KAI) model, is based on the classical statistical th eory 
of nucleation and unrestricted domain growth (lol . 11|. 
For a uniformly polarized FE sample under E ext , the 
KAI model gives the time (i)-dependent change in polar- 
ization AP(i) as 



AP(t) = 2P„[1 - exp{-(t/t ) n }], 



(1) 



where n and to are the effective dimension and character- 
istic switching time for the domain growth, respectively, 
and P s is spontaneous polarization. When the nuclei are 
appearing in time with the same probability, n = 3 for 
bulk samples and n = 2 for thin films In addition, 
to is proportional to the average distance between the 
nuclei, divided by the domain wall speed. Several stud- 
ies have used the KAI model successfully to explain the 
AP(f) behaviors of FE single crystals and epitaxial thin 
films 0. 

Recently, FE thin films have been intensively investi- 
gated for FE random access memory (FeRAM) 0. Most 
commercial FeRAM use polycrystalline Pb(Zr,Ti)0*3 
(poly-PZT) films, and their AP(i) behaviors determine 
the reading and writing speeds of the FeRAM. In such 
non-epitaxial FE films, a domain cannot propagate in- 
definitely due to pinning caused by numerous defects, so 
the KAI model cannot be applied. Therefore, it is impor- 
tant both scientifically and technologically to clarify the 
domain switching kinetics of polycrystalline FE films. 

Numerous studies have examined the AP(t) behaviors 
of polycrystalline FE films, and the reported results vary 
markedly [EM 0,0]. Lohsc et al. measured the polar- 
ization switching currents of poly-PZT films, and showed 
that AP(t) slowed significantly compared to Eq. (1) 
Tagantsev et al. observed similar phenomena for poly- 
PZT films. To explain these behaviors, they developed 



the nucleation-limited-switching (NLS) model. They as- 
sumed that films consist of several areas that have inde- 
pendent switching kinetics: 

AP(t) - 2P s /~Jl - exp{-(t/t ) n }]F(logt )d(logt ), 

(2) 

where P(log to) is the distribution function for log to 
(ij. They assumed a very broad mesa-like function for 
P(log t ), and could explain their AP(t) data. The same 
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FIG. 1: (color online). Schematic diagrams of the pulse trains 
used to measure (a) non-switching polarization (P na ) and (b) 
switching polarization (P 3W ). Time (t)-dependent switched 
polarization AP(t) (c) under various external voltages (Vext) 
at room temperature and (d) under 1.7 V at various temper- 
atures. The dotted and solid lines correspond to fitted results 
using the KAI model and the Lorentzian distribution function 
in log to, respectively. 
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FIG. 2: (color online), (a) Values of n for various T and 
Eext- (b) AP(t) results for (solid symbols) experimental data 
and fitted results using the Lorentzian (solid line), Gaussian 
(dashed line), and delta (dotted line) distributions for log to- 
The inset shows the distribution functions corresponding to 
the fitted results. 



authors also studied La-doped poly-PZT films and found 
that AP(t) at room temperature is limited mainly by 
nucleation, while at a low temperature (T), the switching 
kinetics are governed by domain wall motion, implying 
the validity of the KAI model 

In this Letter, we investigate the polarization switch- 
ing behaviors of poly-PZT films. We can explain the 
measured AP(t) in terms of the Lorentzian distribution 
function for F(log to), irrespective of T. We show that 
such distribution arises from local field variation in a dis- 
ordered system with dipole-dipole interactions. 

Note that (lll)-oriented poly-PZT films with a Ti 
concentration near 0.7 are the most widely used mate- 
rial in FeRAM applications. We prepared our polycrys- 
talline PbZro.3Tio.7O3 thin film on Pt/Ti/Si0 2 /Si sub- 
strates using the sol-gel method. The poly-PZT film 
had a thickness of 150 nm. X-ray diffraction studies 
showed that it has the (lll)-preferred orientation, and 
scanning electron microscopy studies indicated that our 
poly-PZT film consists of grains with a size of about 200 
nm. We deposited Pt top electrodes using sputtering 
with a shadow mask. The areas of the top electrodes 
were about 7.9 xlO -9 m 2 . 

We obtain the AP{t) values of our Pt/PZT/Pt capac- 
itors using pulse measurements @, 0, El, Figure 
1(a) shows the pulse trains used to measure the non- 
switching polarization change (P„ s ). Using pulse Al, we 



poled all the FE domains in one direction. Then, we ap- 
plied pulse A2 with the same polarity, and measure the 
current passing a sensing resistor. By integrating the cur- 
rent, we could obtain the P ns values. Figure 1(b) shows 
the pulse trains used to measure the switching polariza- 
tion (P sw ). Inserting pulse B with the opposite polarity 
between pulses Al and A2, we could reverse some portion 
of the FE domains, so the difference between the values 
of P sw and P ns represents the polarization change due 
to domain switching, namely AP(t). We varied t from 
200 ns to 1 ms, and V ex t from 0.8 to 4 V. The value of 
E ext can be estimated easily by dividing V ex t by the film 
thickness. At T of 80-300 K, we used pulses Al and A2 
with a height of 4 V, which was larger than the coercive 
voltage. Below 80 K, the coercive voltage increases, so 
we increased the pulse height to 6 V [15| . 

Figure 1(c) shows the values of AP(t)/2P s at room 
temperature with numerous values of V ex t- Figure 1(d) 
shows the values of AP(t)/2P s at various T with V ex t 
— 1.7 V. The dotted lines in both figures are the curves 
best fitting Eq. (1). The KAI model predictions deviated 
markedly from the experimental AP(t) values in the late 
switching stage, in agreement with Gruverman et al. [(|. 
In addition, the best fitting results with the KAI model 
gave unreasonable values of n. As shown in Fig. 2(a), 
the values of n varied markedly with T and E ext . In ad- 
dition, in the low E ext region, we obtained n values much 
smaller than 1, which are not proper as an effective di- 
mension of domain growth. Therefore, Eq. (1) fails to 
describe the polarization switching behaviors of our PZT 
films. 

To explain the measured AP(t), we tried simple func- 
tions for F(log to) in Eq. (2). The opposite domain, 
once nucleated, will propagate inside the film, so we fixed 
n=2. The solid circles in Fig. 2(b) show the experimen- 
tal AP(t) at 80 K with V ext = 1.7 V. For P(log t ), 
we tried the delta, Gaussian, and Lorentzian distribu- 
tion functions, as shown in the inset. The dotted line 
indicates the fitting results using Eq. (2) with a delta 
function. Note that this curve corresponds to a fit with 
the KAI model, and thus the classical theory cannot ex- 
plain our experimental data. The dashed line shows the 
Gaussian fitting results. Although this fitting seems rea- 
sonable, some discrepancies occur. The solid black shows 
the fitting results with the Lorentzian distribution: 



F{logto) 



{log to — log ti) 2 -\-w 2 



(3) 



where A is a normalization constant, and w (log £1) is 
the half- width at half-maximum (a central value) [1 €31 ] - 
The Lorentzian fit can account for our observed AP(i) 
behaviors quite well. 

We applied the Lorentzian fit to all of the other exper- 
imental AP(t) data. As shown by the solid lines in Figs. 
1(c) and (d), the Lorentzian fit provides excellent expla- 
nations. Figure 3(a) presents the Lorentzian distribution 
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FIG. 3: (color online), (a) The JSext-dependent Lorentzian 
distribution functions at room temperature, (b) Rescaled 
AP(t) using fitting parameters for the Lorentzian distribu- 
tion function. 



functions used for the 300 K data. As V ex t increases, 
log t\ and w decrease. We rescaled the experimental 
AP(t)/2P s data using (log t - log h)/w. All the data 
merge into a single line, an arctangent function as 
shown in Fig. 3(b). Although not indicated in this fig- 
ure, the experimental data for all other T also merged 
with this line. This scaling behavior suggests that the 
Lorentzian distribution function for log to is intrinsic. 

Note that F(log to) follows not the Gaussian distribu- 
tion, but the Lorentzian distribution. For a statistically 
independent random process, it is a basic statistical rule 
that the resulting distribution should become Gaussian, 
regardless of the process details (l7j . For example, im- 
purities (or crystal defects) inside a real crystal result in 
inhomogeneous broadening of the light absorption line, 
which has a Gaussian line shape. 

However, some studies have observed that magnetic 
resonance line broadening of randomly distributed dipole 
impurities follows the Lorentzian distribution [18j |. The 
first rigorous theoretical result for this problem is that of 
Anderson, who showed that the distribution of any in- 
teraction field component in the system of dilute aligned 
dipoles should be Lorentzian 19|, 2(| • Polycrystalline FE 
films should contain many dipole defects that will act 
as pinning sites for the domain wall motion. To explain 
our observed Lorentzian distribution of log to, we assume 
that a local field E exists at the FE domain pinning sites 



HE) = $ 



E +A 2 



(4) 



where A is the half-width at half-maximum of the E dis- 
tribution function, related to the concentration of pin- 
ning sites. 

In the low E ext region, the domain wall motion should 
be governed by thermal activation process at the pin- 
ning sites. Without E effects, thermal activation results 
in a domain wall speed in the form v oc 1/irj oc exp[- 
(U I k,BT){Eo I E ext )}, where U is the energy barrier and 
Eq is the threshold electric field for pinned domains 21 1 . 
Since E results in a change in the effective electric field 
at pinning sites, the associated to can be expressed as 



to ~ exp 



E { , 



E ext -\-E 



(5) 



Then, the distribution of E results in a distribu- 
tion in log to, using the relation F(logto) = F(E) ■ 
\dE(logto)/d{logto)\. With 



logt x 



UEq 
k B T 



and 



UEq A / 1 



(6) 



(7) 



for 



we can obtain the desired Lorentzian distribution 
F(log t ), i.e., Eq. (3), from Eqs. (4) and (5). 

Our experimental values for log t\ and w agree with 
the analytical forms. Figures 4(a) and (b) plot log t\ vs. 
1/E ext and w vs. 1/E^ xt at various T, respectively. Both 
log ti and w follow the expected £J ea;t -dependence in the 
low E ext region. Note that Eq. (6) is consistent with 
Merz's law 22j, which states that the current coming 
from FE polarization switching should have a character- 
istic time of exp(a/ E ex t), where a is the activation field. 
Using this empirical law, several studies have measured a 
values. For example, So et al. reported a «1700 kV/cm 
for 100-nm-thick epitaxial PZT films Q, and Scott et 
al. re por ted a rj 270 kV/cm for 350-nm-thick poly-PZT 
films [23j. These values are consistent with our room 
temperature value of UE /kBT, i.e., 1400 kV/cm. 

Our model viewed the FE domain switching kinet- 
ics as domain wall motion driven by E ext with a ran- 
dom pinning potential. In the low E ex t region, thermal 
activation at the pinning sites can be important, result- 
ing in the so-called domain wall creep motion. Applying 
atomic force microscopy, Tybell et al. 2l| and Paruch 
et al. 2J] demonstrated that the domain-switching ki- 
netics in epitaxial PZT films are governed by the domain 
wall creep motion. Some theoreticians studied the do- 
main wall creep motion of an elastic string in a random 
potential. They found a linear increase in U with an in- 
crease in T [25(]. The insets in Fig. 4(a) show that the 
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by an external force in presence of a pinning potential 
251 ] . It can be applied to many FE films, since 
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the domain wall motion with a disordered pinning po- 
tential should be the dominant mechanism for AP(t). 
Therefore, the AP(t) studies can be used to investi- 
gate numerous intriguing issues concerning nonlinear sys- 
tems, such as creep motion, avalanche phenomenon, pin- 
ning/depinning transition, and so on. 

In summary, we investigated the polarization switch- 
ing behaviors of (lll)-oriented poly-PZT films and 
found that the characteristic switching time obeyed the 
Lorentzian distribution. We explained this intriguing 
phenomenon by introducing the local electric field due 
to the defect dipolc. 

The authors thank D. Kim for fruitful discussions. 
This study was financially supported by Creative Re- 
search Initiatives (Functionally Integrated Oxide Het- 
erostructure) of MOST/KOSEF. 



FIG. 4: (color online). _E ea; t-dependent (a) log ti and (b) w at 
various T. Note that log t\ and w are proportional to 1/E ext 
and to l/Ee Xt m the low E ex t region, respectively. The insets 
show UEo/ks and A. The solid lines are guidelines for eyes. 

value of UEo/ks obtained from the linear fits in Fig. 4(a) 
increase linearly with T, consistent with the theoretical 
prediction for U [25|. The inset in Fig. 4(b) shows A 
obtained from the fits to Fig. 4(b). Similar exponential 
decay behavior was predicted in a magnetic resonance 
study of randomly distributed dipoles [26| • 

At this point, we wish to compare our model with the 
NLS model. Although both models use Eq. (2), the ori- 
gins and forms for F(log to) are quite different. In the 
NLS model, the FE film consists of numerous areas, each 
with its own and independent to- Subsequently, it was 
suggested that the individually switched regions corre- 
spond to single grains or clusters of grains in which the 
grain boundaries act as frontiers limiting the propaga- 
tion of the switched region [§]. Consequently, the NLS 
model can be applied for polycrystalline films only, and 
the form of F(log to) should depend on their microstruc- 
ture. Conversely, in our model, the interaction between 
dipole defects inside the FE film induces a distribution 
in the local field, which results in F(log to). Therefore, 
both point defects and the grain boundaries could act 
as pinning sites. Using the Lorentzian distribution for 
F(log to), our model can be used for both epitaxial and 
polycrystalline FE films Q. Using Eqs. (2) and (3) with 
small w values, we could successfully explain the AP(i) 
for FE single crystals or epitaxial thin films 0] • We also 
found that our model can explain the AP(i) data for 
poly-PZT films with Ti concentrations of 0.48 and 0.65. 

Note that our model for thermally activated domain 
switching kinetics can be viewed as the famous prob- 
lem that treats the propagation of elastic objects driven 
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